Mon. Not. R. Astron. Soc. 000, 1-7 (2003) 



Printed 2 February 2008 



(MN IMfeX style file v2.2) 



Reconstructing the primordial power spectrum 



S.L. Bridle^*, A.M. Lewis^ J. Weller^ and G. Efstathiou^ 

^Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CBS OH A, UK 
2 CITA, 60 St. George St, Toronto M5S 3H8, ON, Canada 



2 February 2008 



ABSTRACT 

We reconstruct the shape of the primordial power spectrum from the latest cosmic 
microwave background data, including the new results fr'om the Wilkinson Microwave 
Anisotropy Probe (WMAP), and large scale structure data from the two degree field 
galaxy redshift survey (2dFGRS). We tested four parameterizations taking into ac- 
count the uncertainties in four cosmological parameters. First we parameterize the 
initial spectrum by a tilt and a running spectral index, finding marginal evidence for 
a running spectral index only if the first three WMAP multipoles {£ = 2, 3, 4) are in- 
cluded in the analysis. Secondly, to investigate further the low CMB large scale power, 
we modify the conventional power-law spectrum by introducing a scale above which 
there is no power. We find a preferred position of the cut at fee ~ 3 x 10~''Mpc^^ 
although fee = (no cut) is not ruled out. Thirdly we use a model independent param- 
eterization, with 16 bands in wavenumber, and find no obvious sign of deviation from 
a power law spectrum on the scales investigated. Furthermore the values of the other 
cosmological parameters defining the model remain relatively well constrained despite 
the freedom in the shape of the initial power spectrum. Finally we investigate a model 
motivated by double inflation, in which the power spectrum has a break between two 
characteristic wavenumbers. We find that if a break is required to be in the range 
0.01 < fc/Mpc^^ < 0.1 then the ratio of amplitudes across the break is constrained to 
be 1.23±0.14. Our results are consistent with a power law spectrum that is featureless 
and close to scale invariant over the wavenumber range 0.005 fc/Mpc""'^ ^ 0.15, with 
a hint of a decrease in power on the largest scales. 

Key words: cosmology:observations - cosmology:theory - cosmic microwave back- 
ground - large scale structure 



1 INTRODUCTION 

Measurements of the Cosmic Microwave Background (CMB) 
radiation have taken a leap forward with the recent an- 
nouncement of the findings of the Wilkinson Microwave 
Anisotropy Probe (WMAP). With these new data it is pos- 
sible to set important constraints on the shape of the pri- 
mordial power spectrum and hence to begin to differentiate 
between the plethora of models for the early universe. One of 
the most intriguing results comes from a a combined anal- 
ysis of the WMAP data with the two degree field galaxy 
redshift survey (2dFGRS) data (Spergel et al. 2003; Peiris 
et al. 2003), which indicates that the primordial power spec- 
trum might have curvature. The addition of Lyman-a forest 
data on smaller scales to strengthen this conclusion is how- 
ever contentious (Seljak et al. 2003). The low quadrupole 
and octopole observed in the CMB temperature power spec- 
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trum (Spergel et al. 2003) has a low probability in standard 
models, and may be an indication of some feature in the 
initial power spectrum on very large scales. 

Both inflationary Big-Bang (Guth 1981; Linde 1982; 
Albrecht and Steinhardt 1982) and more speculative cyclic 
ekpyrotic (Steinhardt and Turok 2002) models of the early 
universe predict very nearly Gaussian scalar perturbations 
in the primordial radiation dominated era. The shape of the 
perturbation power spectrum depends on the exact model, 
which typically involves various unknown parameters. The 
objective of this paper is to constrain the shape of the ini- 
tial power spectrum directly from observational data with 
as few assumptions as possible. 

A wealth of cosmological information can be obtained 
from analysing the shape of the cosmic microwave back- 
ground (CMB) radiation temperature fluctuation power 
spectrum. However it has been shown that the efi'ect of 
changing the cosmological parameters can be exactly mim- 
icked by changes in the shape of the primordial power spec- 
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trum (Souradccp ct al. 1998; Kinney 2001). By including 
measurements of the CMB polarization and the late time 
matter power spectrum, as measured for example by galaxy 
redshift surveys, the degeneracy cam be broken because the 
cosmological parameters affect these data in different ways. 
In this paper we combine temperature and polarization data 
from the WMAP observations and other CMB observations 
on smaller scales with constraints on the matter power spec- 
trum from the 2dFGRS (Percival et al. 2002). 

Since the initial power spectrum is an unknown func- 
tion one is forced to parameterize it. There are numerous 
possibilities. Wang et al. (1999), Wang and Mathews (2002) 
and Mukherjee and Wang (2003a) use a number of bands 
in wavenumber. Mukherjee and Wang (2003b, c, a) also use 
a model independent approach but using wavelets. Barriga 
et al. (2001) test a particular inflationary model, in which a 
phase transition briefly halts the slow roll of the inflaton. 

Inflationary models generically predict a monotonically 
slowly varying power spectrum, determined by the shape 
of the inflationary potential. More complicated inflationary 
models can give more interesting spectra at the expense of 
introducing parameters that are fine tuned to give effects in 
the small range of observable wavenumbers. In this paper 
we consider various different power spectrum parameteriza- 
tions, motivated by theoretical models or features of the ob- 
served data. We also adopt a model independent approach 
which allows general trends or unexpected features to be 
detected. 

In Section 2 we specify the framework used in the rest 
of the paper. We summarize the results of the conventional 
power spectrum parameterization (ie. a power law slope and 
running spectral index) in Section 3, and discuss the possi- 
ble implications of the small large scale power detected by 
WMAP in Section 4. Section 5 describes a reconstruction 
of the primordial power spectrum in wavenumber bands on 
sub-horizon scales. We test a double inflation model over a 
particular range of wavenumbers in Section 6. In each sec- 
tion we discuss the implications for estimates of the other 
cosmological parameters. 



2 FRAMEWORK 

The primordial scalar power spectrum is defined by 

(IxP) = / dlnfc'P^(fc), where x is the super-horizon comov- 
ing curvature perturbation in the early radiation dominated 
era. The commonly assumed power-law power spectrum pa- 
rameterization is then 

Here ns{k) = d\n'Pxik)/d\nk +1 is the conventional def- 
inition of the scalar spectral index, where rig = 1 corre- 
sponds to a scale invariant power spectrum (we use feso = 
0.05Mpc~^ throughout). The power spectrum amplitude Ag 

1/2 

determines the variance of the fluctuations, with Ag ~ 
5 X 10~^ to give the observed CMB anisotropy amplitude. 

In slow roll inflationary models, it is expected that the 
spectrum varies only very slowly and that |ns — 1| is much 
smaller than unity (Lyth and Riotto 1999). In general there 
is a direct relation between the potential of the inflaton field 
and the spectral index. As the potential evolves during infla- 



tion the spectral index can vary slightly as different modes 
leave the horizon. This can be characterized by including a 
second order term in the logarithmic expansion of the power 
spectrum rirun, defined by 

ln7'.(.)=ln^. + (n.-l)ln(A)+!^(,,(_L))\(2) 

The value of Us therefore depends on the pivot scale used, 

for example to convert to a new pivot scale the relation is 
f^sik'so) = 7is(fcso)+nrunlog(fcso/fcso)- More generally double 
inflation models or multiple field inflation can lead to breaks 
and spikes in the primordial power spectrum, see eg. Linde 
(1990). This motivates a more general parameterization of 
the primordial power spectrum in terms of amplitudes over 
discrete bands in wavenumber. 

The primordial power spectrum is related to the linear 
CMB anisotropy power spectrum Cf^ by a transfer func- 
tion Tl^{k) via 

Cf" oc j din kVx{k)Te''{k)Tlik). (3) 

where X and Y denote the various temperature and polar- 
ization power spectra. The transfer function for a mode of 
wavenumber k peaks at a multipole of about £ ~ kdA where 
dA is the angular diameter distance. For the concordance 
model the relation is roughly £ ~ (le4Mpc)fe. However the 
detailed shape of the transfer function is complicated and 
a range of wavenumbers contribute to each multipole (see 
eg. (Tegmark and Zaldarriaga 2002)). 

In this paper we vary four cosmological parameters, 
using fiat priors on the baryon density fib/i^, the cold 
dark matter density Q,ch^, the Hubble constant h — 
//o/( 100km s~^Mpc~^), and the redshift of reionization ZeS 
(we assume 6 < Zgs < 30). The true ionization history may 
be complicated, so we assume there exists an effective red- 
shift 2eff at which a single rapid reionization gives a good 
approximation to the true CMB anisotropy. The tempera- 
ture CMB anisotropy is fairly insensitive to the details of 
reionization, and the effect on the polarization is only on 
large scales where it is largely hidden by cosmic variance and 
the current observational noise, so this assumption should 
not affect our results significantly. Throughout we assume 
that the universe is fiat with a cosmological constant. We 
assume purely Gaussian adiabatic scalar perturbations and 
ignore tensor modes in this paper. Since tensor mode per- 
turbations decay on sub-horizon scales they only affect the 
large scale CMB anisotropy. Given our assumptions, if we 
find excess power on large scales this could, equivalently, be 
due to a tensor contribution rather than the shape of the 
scalar initial power spectrum. 

We use the latest WMAP^ (Verde et al. 2003; Hinshaw 
et al. 2003; Kogut et al. 2003) temperature and temperature- 
polarization cross-correlation anisotropy data. We also in- 
clude almost independent bandpowers on smaller linear 
scales (800 < £ < 2000) from ACBAR^ (Kuo et al. 2003), 
CBI (Pearson et al. 2002) and the VSA (Grainge et al. 2002). 

We constrain the matter power spectrum at low red- 
shift by using the galaxy power spectrum measurements of 

^ http : //lambda.gsfc.nasa.gov/ 

^ http: //cosmology.berkeley.edu/group/swlh/acbar/data 
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Figure 1. Marginalized distributions of the running spectral in- 
dex slope parameters including (solid) and without (dashed) the 
WMAP temperature data at £ < 5. 

the 2dFGRS (Percival et al. 2002) over the range 0.02 < 
k/{hMpc~^) < 0.15. We assume that the galaxy power spec- 
trum measured by the 2dFGRS is a simple unknown multi- 
ple of the underlying matter power spectrum (linear bias), 
so in our analysis the 2dFGRS data serves to constrain the 
shape but not (directly) the amplitude of the matter power 
spectrum. In principle the matter power spectrum is directly 
proportional to the primordial power spectrum however in 
practice, due to the finite volume observed, the data con- 
strain a smeared version of the underlying matter power 
spectrum. This makes it harder to detect any sharp features 
in the primordial power spectrum, as investigated by El- 
garoy et al. (2002). 

To evaluate the posterior distributions of parameters 
given the data we use the Markov-Chain Monte Carlo 
method to generate a list of samples (coordinates in pa- 
rameter space) such that the number density of samples is 
proportional to the probability density. We use a modified 
version of the CosmoMC^ code, using CAMB (Lewis et al. 
2000) (based on CMBFAST (Seljak and Zaldarriaga 1996)) 
to generate the CMB and matter power spectrum transfer 
functions. CosmoMC uses the Metropolis-Hasting algorithm 
to explore the posterior probability distribution in a piece- 
wise manner, allowing us to exploit the fact that for each 
transfer function computed many different values of the ini- 
tial power spectra parameters can be changed at almost no 
additional computational cost. For further details see Lewis 
and Bridle (2002); Christensen et al. (2001); Verde et al. 
(2003) and references therein. Most of the chains for the 
analysis here were computed in around a day using spare 
nodes of the CITA beowulf cluster, with each node running 
one chain parallelized over the two processors. Between 4 
and 20 converged chains were generated for each set of pa- 
rameters, burn in samples were removed, leaving of the order 
of 10^-10® accepted positions in parameter space from which 
the results in this paper were computed. 



3 POWER LAW SPECTRA WITH AND 

WITHOUT A RUNNING SPECTRAL INDEX 

In most parameter studies the scalar initial power spectrum 
is parameterized by a constant spectral index Us and an am- 
plitude. Even with only two parameters defining the primor- 
dial power spectrum there are large degeneracies between 
these and the other cosmological parameters. For example 

^ http : //cosmologist.info/cosmomc/ 



larger baryon densities decrease the height of the CMB sec- 
ond acoustic peak thereby mimicking the effect of a red spec- 
tral tilt. 

Using the WMAP (Hinshaw et al. 2003; Kogut et al. 
2003) power spectrum results alone we find a tight marginal- 
ized constraint - SOAiflth^ ~ 0.025) = 1.04 ± 0.02 (68% 
confidence). However the orthogonal direction is very poorly 
constrained with" -I- 30. 4(^5/1^ - 0.025) = 1 ± 0.14. The 
best fit model to the WMAP data (Spergel et al. 2003) with 
Us = 0.97 can be closely approximated by completely dif- 
ferent models with Us > 1.1, higher reionization redshifts, 
rapid Hubble expansion and high power spectrum ampli- 
tude. On integrating out the value of Qbh^ the constraint 
on the spectral index weakens to Us = 1.05 ± 0.08. Sim- 
ilarly the amplitude of the primordial curvature fluctua- 
tions on 0.05Mpc~^ scales is constrained by WMAP to be 
AI^^ = (5.5 ± 0.8) X 10"^, whereas the parameter combi- 
nation A^^e"^ = (4.3 ± 0.1) X lO"'^ is much more tightly 
constrained since it comes directly from the observed tem- 
perature anisotropy amplitude. To partially break these de- 
generacies the WMAP team adopt a prior on the reioniza- 
tion optical depth of r < 0.3. By adding the 2dFGRS and 
additional CMB data at £ > 800 we find the parameter con- 
straints which broadly agree with the analysis reported by 
the WMAP team (Spergel et al. 2003), see Fig. 5 below. 

By adding a running spectral index we find the 
marginalized 1-sigma result n^un = —0.04 ± 0.03 shown by 
the solid line in the right hand panel of Fig. 1, in rough 
agreement with the WMAP team. For the pivot scale used 
(feso = 0.05Mpc~'^) a red tilt is preferred (left hand panel of 
Fig. 1). As expected, the effect of adding nrun as a free pa- 
rameter is to increase the uncertainties on the cosmological 
parameters but only within the original uncertainties. 

We flnd that the evidence for running comes predomi- 
nantly from the very largest scale multipoles. When we ex- 
clude e < 5 multipoles from the WMAP temperature (TT) 
likelihood the running spectral index distribution shifts to 
becomes highly consistent with nrun = 0, as shown in Fig. 1 
. The constraints on all cosmological parameters except for 
ris and nmn are changed very little on removing the lowest 
multipoles. The running parameterization is therefore not 
ideally suited to the data: WMAP provides some evidence 
for low power on the very largest scales, but this is only 
crudely fit by using a running spectral index. The WMAP 
analysis relies on Lyman-a forest at wavenumbers greater 
than k ~ O.lMpc"^ to give evidence for more red tilt on 
small scales consistent with a running index, but the valid- 
ity of this analysis is in serious doubt Seljak et al. (2003). We 
conclude that the marginal preference for a running spectral 
index from our CMB -f- 2dFGRS analysis is primarily driven 
by the first three CMB multipoles. 



4 POWER SPECTRUM ON THE LARGEST 
OBSERVABLE SCALES 

The quadrupole and octopole estimators observed by 
WMAP are low compared to the other large scale multi- 

* This wide spread may be partly due the to approximations used 
in the WMAP likelihood parameterization. 
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Figure 2. Top: Marginalized probability distribution of a large 
scale power cut-off parameter kc for which Vxik < kc) = 0. Lower 
panels: a concordance model with k^ = (solid lines) and with 
kc = 2.7 X 10"* (dashed lines) compared to the WMAP data. 
Dark and light error bars show the 68 and 95 confidence limits on 
the theoretical value at each multipole (for the TE error bars we 
use the Gaussian contribution to the likelihood given the observed 
TT estimator, taking into account the correlation between the TT 
and TE estimators, assuming the WMAP best fit model for Cj"^ 
and Cf^, where (Cj^f < Cj^^Cf^). 



poles. In a given model the low multipole estimators have 
a wide x^-like distribution (which has the peak below the 
ensemble mean), so the low values could just be chance. 
However any model that predicts small values for the low 
multipoles would be favoured by the data, by a factor of up 
to about fourty, so this could be a hint that there is a sharp 
fall in power on the largest scales. Tegmark et al. (2003) find 
that the quadrupole and octopole appear to be aligned, per- 
haps indicative of some anisotropic effect on on very large 
spatial scales which would not be well modeled by a sta- 
tistically isotropic power spectrum model. Here we assume 
that the alignment is a coincidence, and proceed to consider 
whether the shape of the initial power spectrum could help 
to explain the low large scale signal. 

There is only a limited amount of information on the 
largest scales due to cosmic variance, so we cannot hope 
to fit many extra parameters. We choose to assume there 
is a sharp total cut-off in the primordial power on scales 



larger than k = kc (for a discussion of an exponential cut- 
off in closed models see Efstathiou (2003); see Contaldi et al. 
(2003) for possible motivation for a cut-off). As discussed in 
the previous section a constant spectral index is a good fit on 
smaller scales, so we parameterize the primordial spectrum 
as 







(4) 



Marginalizing over the other parameters we find the con- 
straint on kc shown in the upper panel of Fig. 2. We find 
a preference for a cut-off at kc = {2.7t1'l) ^ 10~*Mpc~^, a 
scale which can give a significantly lower quadrupole and oc- 
topole without significantly affecting the higher multipoles 
(lower panels of Fig. 2). However, according to this model 
a spectrum with no cut-off is not strongly excluded by the 
data. Our parameterization cannot achieve values for the 
quadrupole as low as observed because there is a signifi- 
cant Integrated-Sachs Wolfe contribution to the quadrupole 
from scales smaller than the cut. The best-fit cut-off model 
has a very similar probability to the best-fit running model, 
although the cut-off model is marginally preferred. The con- 
straints on the cosmological parameters are virtually unaf- 
fected by adding this free cut-off scale. 

A comparison with Spergel et al. (2003) 's assessment of 
the significance of the low CMB power at small multipoles 
is not straightforward because the statistical tests differ. A 
more detailed analysis of this important issue is required, 
but given the results of Fig. 2 it would seem premature to 
discard simple continuous power law models. 



5 POWER SPECTRUM RECONSTRUCTION 
IN WAVENUMBER BANDS 

On smaller scales there is now a wealth of data available 
to constrain the primordial power spectrum in considerable 
generality. Here we choose to parameterize the spectrum in a 
number of bands, which has the advantage that it is free from 
any assumptions about the underlying model. This general- 
ity can in principle reveal unexpected features in the pri- 
mordial power spectrum on scales comparable to (or larger 
than) the band width. The disadvantage is that there are a 
large number of free parameters so care is needed to obtain 
meaningful results. 

The initial power spectrum is only indirectly con- 
strained by the observational data. Since we do not model 
the biasing in the 2dFGRS galaxy power spectrum, the only 
direct constraint on the amplitude (as opposed to the shape) 
comes from the CMB anisotropy power spectra. On scales 
smaller than the horizon size at reionization {£ ~ 40) the 
observed power in the CMB anisotropy scales as e~^^'P^(fe), 
depending on the optical depth to reionization r. We there- 
fore parameterize the shape of the initial power spectrum by 
the value of e~^'^Vx{k) at a set of points, and linearly inter- 
polate (in k) between the points. We then refer to the power 
at each point as the power in a 'band' at that wavenumber. 
Our parameters bi are defined by 



(fei+i-fe)()i+(fc-fei)i)i 



b„ 



h <k < ki+i 



(5) 



where the first line applies for 1 i ^ n — 1. We set the 
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Figure 3. CMB power spectra for top hat primordial power spec- 
tra at the positions of the bins used for Fig. 4. The top line is the 
sum of all the lower lines, ie. the usual CMB power spectrum. 



position of the first band to fci = and that of the second 
band to k2 = 0.005Mpc~^, so we only have two bands over 
the region where the data are effected by the modes which 
are super-horizon at reionization and hence do not scale sim- 
ply with e"^^. Subsequent band positions on smaller scales 
are logarithmically spaced with ki+i = 1.275 ki, where the 
constant is chosen such that fcie — 0.15Mpc~^. We assume a 
flat prior on each of the bi . This prior is equivalent to a flat 
prior on the underlying power spectrum e^^bi, together with 
a prior on the optical depth P{t) oc e~^"^. Constraints on 
quantities sensitive to r depend strongly on the choice prior 
for large numbers of bands. This is one reason we reconstruct 
bi , which are more directly constrained by the data and have 
the main dependence on r taken out, rather than Vxik) di- 
rectly. Our choice of prior gives a posterior constraint on the 
optical depth similar to (but somewhat broader than) that 
with the simple spectral index parameterization. 

Fig. 3 shows the CMB temperature power spectra for 
top hat primordial power spectra which are zero except for 
the range (ki-i + ki)/2 < k < {ki + fci+i)/2 (except for 
i — 1 where < fc < ^2/2) for a concordance model. All the 
primordial power spectrum top hats have the same absolute 
height. The result illustrates the rough correspondence be- 
tween the position of the bands we use and the CMB power 
spectrum peaks, showing that there is one band covering 
much of the large scales, and several bands over the first 
acoustic peak, with only one or two bands for each of the 
subsequent few peaks. 

Our reconstruction of the shape of the initial power 
spectrum is shown in Fig. 4. The crosses and error bars show 
the peak and 68% upper and lower limits of the marginalized 
distributions. The probability distribution of all but the last 
band is close to Gaussian, although some bands are corre- 
lated at around the fifty per cent level. The maximum corre- 
lation is 80 per cent, which is a positive correlation between 
bands around the first acoustic peak of the CMB. We have 
checked that our constraints on the sub-horizon bi agree well 
with those obtained by fixing the optical depth to r = 0.17 
rather than marginalizing over it. The overall shape is con- 
sistent with a featureless scale invariant power spectrum, 
though perhaps an overall red tilt is discernible, bearing in 
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Figure 5. Constraints on cosmological parameter estimates after 
marginalising over the power spectrum amplitudes in 16 bands 
(solid lines), and for comparison equivalent results in the constant 
spectral index model (dashed lines). Contours enclose 68% and 
95% of the probability. 

mind the large error bar on the small scale point. We do not 
see the feature at A; ~ 0.015Mpc~^ suggested by Mukherjee 
and Wang (2003b). Our reconstruction is complementary to 
that in Mukherjee and Wang (2003a) since our inter-band 
separation is one third of theirs. Also, we use linear interpo- 
lation instead of wavelets or top-hat bins, which makes the 
primordial power spectrum relatively smooth. 

In Fig. 5 we investigate how robust the cosmological 
parameter estimates are to this general primordial power 
spectrum. Dashed lines show the cosmological parameter 
constraints assuming a power law primordial power spec- 
trum, marginalising over n^. Overlaid, the solid lines show 
the result after allowing freedom in the amplitudes of the 
16 bands. The error contours are broader but remarkably 
show that it is still possible to recover interesting constraints 
on the cosmological parameters even if one allows the large 
amount of freedom in the primordial power spectrum shape. 

It is important to emphasize that this method of power 
spectrum differs from earlier work such as that by Gawiser 
and Silk (1998) because the MCMC reconstruction prop- 
erly accounts for the allowed ranges of all of the parameters 
defining the model. 



6 BROKEN POWER SPECTRUM 

In this section we use a parameterization motivated by dou- 
ble infiation, previously explored by Barriga et al. (2001). 
In this model 

{A k < ki 

Ck°"^ fci < fc < fca (6) 
B k> k2 . 

where C and a are chosen such that the power spectrum 
is continuous. As in Barriga et al. (2001) we take flat pri- 
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Figure 4. Reconstruction of the shape of the primordial power spectrum in 16 bands after marginalising over the Hubble constant, 
baryon and dark matter densities, and the redshift of reionization. 
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Figure 6. Marginalized constraints on the parameters of the bro- 
ken power spectrum model and the correlations with the baryon 
density. The constraint on t is changed relatively little on allow- 
ing the 16 bands, becoming r = 0.13 ± 0.07 
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Figure 7. Primordial power spectra of random samples from the 
broken power spectra chains. 



ors on fci, A/B and log(fc2/fci) and we limit the parame- 
ter space to 0.01 < fci/Mpc"^ < 0.1, 0.3 < A/B < 7.2 
and 0.01 < log(fc2/fci) < 4. In addition we impose a prior 
fc2 < O.lMpc^^ so that we explore a transition only in the 
region probed by the observational data. 

In Fig. 6 we show the constraints on the model pa- 
rameters after marginahsing over four cosmological param- 
eters. A conventional scale-invariant spectrum corresponds 
to A/B = 1 and we can see that this possibility is very 
close to the la contours. Values of A/B higher than unity 
are preferred, corresponding to a drop in the initial power 
spectrum on going from large to small scales. This is not 
surprising since the data favour red tilts, and in this pa- 
rameterization a tilt is obtained by having a long transition 
between a higher flat spectrum on large scales and a lower 
flat spectrum on small scales. The distribution is slightly bi- 
modal, with the first mode roughly corresponding to a wide 
transition straddling the first CMB acoustic peak and the 
second corresponding to a drop in power in the dip between 
the first and second peaks. The expected strong correlation 
between A/B and the baryon density is clear. We find that 
removing the first three multipoles has no effect on the pa- 
rameter constraints on this model. 

We show in Fig. 7 the power spectra corresponding to 20 
samples in the MCMC chain. There is considerable spread, 
but the transition, if any, occurs at around k ^ 0.04Mpc~'^. 
This corresponds roughly to the scale at which the 2dF- 
GRS data becomes statistically significant. Since the 2dF- 
GRS data may disfavour a transition at higher wavenumber 
and the WMAP data would disfavour a transition at lower 
wavenumber then neither the amplitude or the scale of the 
break is likely to be of any physical significance. 

The effect of adding these additional parameters on the 
estimates of cosmological parameters is small; the biggest ef- 
fect is to widen the allowed range of Qhh^ to include smaller 
values. This is due to the fact that a preferred position for 
any transition is between the first two acoustic CMB peaks, 
and it is also this ratio that determines the baryon density. 



7 CONCLUSIONS 

We have explored various parameterizations of the primor- 
dial power spectrum and found that in each case a simple 
scale invariant spectrum is an acceptable fit to the data. 



© 2003 RAS, MNRAS 000, 1-7 



Reconstructing the primordial power spectrum 7 



Some deviation towards a red tilt is preferred and there is 
a marginal indication of a cut-off in power on the largest 
observable scales. Significantly, even if the power spectrum 
is allowed to vary over a large number of wavenumber bands 
on sub-horizon scales the reconstructed spectrum is found 
to be featureless and close to scale invariant. In addition we 
have shown how the constraints on the other cosmological 
parameters are affected by the addition of this extra freedom 
in the primordial power spectrum shape: the error bars on 
cosmological parameters are significantly amplified but the 
resulting constraints are still strong and roughly similar to 
those before the recent WMAP results. The unexpectedly 
small power on superhorizon scales observed by WMAP, if 
real, provides marginal evidence for a drop in the primordial 
power spectrum on the largest observable scales. 
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